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Abstract 

We investigate the thermal behaviour of gluon and chiral condensates within an effective 
Lagrangian of pseudoscalar mesons coupled to a scalar glueball. This Lagrangian mimics 
the scale and chiral symmetries of QCD. 

1. INTRODUCTION 

The construction of effective models for low energy QCD is strongly constrained by global 
symmetry aspects. The most important one is the chiral symmetry which, together with 
its spontaneous breakdown at low energies, was therefore already extensively studied. 
Another symmetry of the classical QCD-Lagrangian in the limit of vanishing quark current 
masses is the dilatation symmetry or scale invariance. This symmetry exhibits an anomaly, 
i.e. it is broken on the quantum level by radiative corrections [Q. 

An effective realisation of the QCD scale anomaly was found in the early 80's f^, by 
adding to the classical Lagrangian a scalar glueball field X with an interaction potential 
of the form 

V{X) = CX^ [In {XjXp) - 1/4] . (1) 

Here C and X^, are parameters to be specified later. This leads to = —CX^ for the trace 
of the energy-momentum tensor which is also the divergence of the dilatation current [^ . 
This quantity with scaling dimension 4 can then be identified with the trace anomaly in 
QCD, 

^^,QCD^ = {^G;^G^n ^ -{CX') . (2) 

A revival of this idea came with the articles of Campbell et. al. @] who, somewhat 
controversially, suggested to regard {X) as an order parameter for a deconfining phase 
transition. Similar to the restoration of spontaneous broken symmetries at high tempera- 
ture and/or densities, deconfinement should signal itself in the effective potential (EP) of 
the A'-field by a phase transition from {X) = A'o 7^ to {X)t^ = 0. 

The breaking of the two symmetries in low-energy QCD mentioned above is connected 
with the appearance of two important vacuum condensates, the chiral condensate (qq) and 
the gluon condensate {G^^G"'^'') , respectively. Therefore these expectation values have 
been considered as order parameters for symmetry restoration. 
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Campbell et al. conjectured a strong correlation between the behaviour of these two 
condensates at rising temperature or density. The main relation, under the assumption of 
factorisation, is 

(99) «(^) (Tr(f/ + f/t)), (3) 

with U denoting the chiral matrix field, defined below. This means that a vanishing of 
{X), i.e. of the gluon condensate, drives also the chiral condensate to zero. Therefore the 
important question is whether the temperature scale for vanishing of the chiral condensate 
is dominated by the intrinsic chiral dynamics, which drives the disappearance of (Tr(f/ + 
t/^)), or by the vanishing of {X). 

The idea in 0] has triggered a lot of (sometimes rather phenomenological) work inves- 
tigating the scaling behaviour of low energy Lagrangians at finite temperature or density 
0]-|T^. Here, we will follow the original suggestion in [Q to determine the thermal be- 
haviour of {X) via the minimum of the effective potential, by performing a systematic 
calculation of the latter. We will see that the conjectures in about the character of the 
chiral and deconfinement phase transitions are only partially supported by our detailed 
calculation. 

We start from the nonlinear cr-model, extended by a scalar glueball field in order to 
mimic the scaling properties of QCD |Q. For this effective low-energy Lagrangian we 
calculate in Sec. 2 the Gaussian effective potential (GEP) for the dilaton field X. From 
the GEP we can read off the temperature dependence of the expectation value {X), leading 
via the identification (|^) to the thermal behaviour of the gluon condensate. This will be 
done in Sec. 3, where we also discuss the influence of thermal excitations of the A'-field 
on the chiral condensate. In Sec. 4 we will discuss our results. 



2. THE GAUSSIAN EFFECTIVE POTENTIAL 

We start with the same Lagrangian as in [Q, but neglecting the contribution of the rj' 
meson: 
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= ZiL _ Tr \dM^d>'U 
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MgiU + U^)]+^d^Xd^'X-V{X). (4) 



f/ = exp i ^ A" 



a=l 



Mq = diag (m„, m^, m^, . . .) , Xq = {X)t= 



(5) 



A" are the generators of flavour SU(A^) (A^ is the number of quark flavours), and V{X) 
was defined in (|l]). 

For X = const. = Xq this Lagrangian reduces, up to an additative constant, to the non- 



linear (j-model [|I3|. As long as the explicit chiral symmetry breaking by the quark mass 
term A4q is small, terms containing higher derivatives of U can be neglected in comparison 
with the first term in M), and (HI) is a good approximation to the chiral dynamics of the 
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A^^ — 1 Goldstone bosons of QCD at low energy. We will thus consider the cases N = 2 
and N = 3. 

We want to calculate the EP Ves{'^) with X = {X). We expect a breakdown of the 
model in the limit X 0: The expansion of the logarithmic tree level potential V{X) 
(0) in a power series near X = leads to singular coefficients for powers larger than 3. 
At the origin we will also pick up singularities in coupling constants related to the chiral 
field U, because of the effective scaling of all ?7-terms in with powers of X (see below). 
This means that the theory cannot be expanded into a Taylor series around X = 0. As 
mentioned in [^] a particle interpretation at the origin is thus not possible. Campbell et 
al. interpreted this phenomenon as indication that in the phase where X = the physical 
degrees of freedom of the model, the Goldstone bosons and the glueballs X, are no 
longer relevant, and that this phase corresponds to a deconfined phase where these mesons 
have dissolved into their quark content. 

The EP implements radiative corrections to the tree level potential V{X) arising from 
virtual and (at non-zero temperature) real excitations of the degrees of freedom in the 
Lagrangian. Since the latter are no longer well-defined in a state where X = 0, the EP 
cannot be perturbatively calculated near X = 0. Thus a direct comparison of the effective 
potential (= free energy) at the origin X = and at the local minimum of Ves{X) is not 
possible in any perturbative approach. Therefore, it is for example impossible to decide 
whether the phase transition connected with the disappearance of the scalar condensate 
{X) is of first or second order. This point was apparently missed in Ref. [Q. However, 
we can still calculate the shift of the local minimum of l4ff(A?) at A" 7^ as a function 
of temperature and investigate at which temperature scale visible thermal effects set in 
which could indicate the disappearance of the gluon condensate. 

We also encounter another familiar problem connected with the EP in the case of 
broken symmetries. The standard loop expansion |]14| for the EP leads to complex contri- 
butions for values of X, for which the tree level potential V{X) is concave, i.e. its second 
derivative is negative. This problem has triggered a large amount of work, ranging from 
an interpretation of the imaginary part as a sign for instability for the corresponding X 



state [|T5[ to techniques which avoid the complex contributions altogether [^, |T] 

In |T3[ values of X for which Ves{X) was complex were interpreted as false vacua, and 
the imaginary part was related to the decay rate per unit volume of the false vacuum. 
One knows, however, from the definition of the effective potential that the exact EP is 



real [|1^]. Therefore we regard the appearence of complex contributions as a sign for an 
unsuitable expansion scheme. We are thus led to choose an approximation where all 
expansion coefficients are real. 

We choose a variational ansatz, the Gaussian effective potential (GEP), which was ex- 



tensively studied by Stevenson et al. |jT9[ . Most other methods, which avoid the imaginary 
part, are only valid in distinct temperature regions. The interpolated loop expansion 
is only tractable for low temperatures while resummation techniques are valid only for 
high temperatures [1^. With the GEP ansatz we are able to cover the whole temperature 
region. A more detailed discussion of the problem of the imaginary part and a comparison 
of different methods to solve the problem is given in 
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The principle of the GEP is the introduction of an arbitrary mass parameter f2 for 
the dilaton field X. The effective potential is then be minimized with respect to this 
parameter, allowing only for real values of Q. This is in analogy to the Ritz variational 
principle in quantum mechanics, where the ground state energy is the minimal energy 
under variation of the wave function. We will introduce the variational parameter Q later 
within the path integral formalism by a method developed by Okopihska ||2^ . 

We start from the finite temperature partition function Z^. After expanding the chiral 
field U, Eq. (|^), in terms up to fourth order in the Goldstone fields $ or 9^$, we get for 
the partition function 



a=l 



Zf^ir, K] = Af n ^[-^i exp { i / d^x £"^(<l>^ 9/,<l'^ X, d^X) + J°<l>° + KX 



A/" is a normalisation constant, and C^^ is given as 
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(7) 

J° and ii' are external sources for the Goldstone bosons and dilaton fields, respectively. 
The first term is independent of Therefore we add this term to the potential V{X) 
and define 



V^{X):=V{X)-g;,X-' with g;, -.= l^Y.^^ ■ 

^0 q 



We will use the real-time formalism. Therefore the integration of Xq is along the real- 
time path C in the complex plane [^. In order to evaluate we will use the sattle point 
approximation. This means that we expand the dilaton field X around the solution of the 
classical equation of motion 



dC{x,u^) dC{x,u^) 

" did.X^i) dX^i 



(9) 



4 



In we substitute the field $ by its expectation value ($(a;)) = 0. Further we define the 
dilaton fiuctuating field X := X — X^i and transform the measure for X to V[X] = T>[X]. 
Because of the starting Lagrangian @) which contains no X independent terms the kinetic 
term for the meson fields scales with [Xd/ X^f". 



<I>kin 



2 \X^ 



Xn 



(10) 



where is the meson mass in the non-linear sigma model (see eq. (|26D below). In order 
to obtain a standard kinetic term for we transform the meson fields by 



11) 



This transformation leads to an infinite Jakobi determinant for the integration measure of 
the $ field, 



Vl^"] — > X Det 



:i2) 



which we absorb in the normalisation J\f. After these manipulations we can write the 
partition function as 



z^[r,K]=Mj X{v[^'']v[x] 

a=l 
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^0 afecd 

At this stage of the calculation we have for the mass of the A'-field 



(19) 



(20) 



The meson masses scale with ^^JX^i/Xq, and we have defined 
Further we combine the interaction terms by 

-^int = -^^int + ^Xint + '^A'Sint • 



(21) 



(22) 



We now introduce the variational parameter f2 by the method developed in This 
is done by adding and subtracting the quadratic term Vl? X"^ /2 to the Lagrangian (|1^). We 
now regard — X"^ /2 as a new mass term for the A'-field and get therefore an additional 
contribution to the interaction part 



C 



int 



c,,,,-]^{Ml-n'')x\ 



(23) 
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This corresponds to a two point vertex, indicated by a cross in the Feynman diagrams (see 
Fig. 1). We thus have to evaluate the expression 



Z^[J\K]=Uf \{V[^'']V[X] 

a=l 



(24) 



exp 




This expression (|2^) is by construction independent of the arbitrary parameter fi. But 
in practice, in order to evaluate (p^), we have to expand the exponential containing the 



interaction part Ci^^ and truncate the series at some point. This truncation introduces a 
dependence on Q. Given such an approximation to the EP as a function of Q, one then 
uses the principle of minimal sensitivity with respect to the arbitrary parameter Q to fix 
^1{X) as a function of X by the condition 



. 



(25) 



In our present work we expand exp{i / d^xC^^} up to first order. This means that we 
only include vacuum bubbles containing at most one vertex. The Feynman diagrams with 
one vertex which appear in our calculation are shown in Fig. 1. 

We investigate three cases. Case 1 corresponds to a pure glueball scenario, neglecting 
in (P all ?7-field terms. In case 2 we include a SU(A^) flavour representation with N equal 
masses for the [/-fields. In case 3 we consider a broken SU(3) flavour representation with 
fTT's 7^ ^n,d, dcuotcd by SU(3)b. For the last case the physical meson masses in zeroth 
order are 
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c{mq + nis) = ml = ml = m^ = m^ 



m^= 2c{mq + 2ms) /?) = 



(26) 



where = m„ = m^. 

The result for the partition functions for the three cases, neglecting the normalisation, 

is 



= exp < 
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(29) 



The integrals /„ are defined in the appendix. From the functional Z[K] = exp(iiy[ii']) we 
get the EP by first making a Legendre transformation 



with 



T\X{x)] = W[K] - j d!^xK{x)X{x) , 
6W[K] 



X{x) 



5K{x) 



(30) 
(31) 
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The expansion of the effective action r[A:'(x)] in powers of derivatives 



XI 



(32) 



contains as its first term the EP. 

We apply this procedure to the partition function (P7|)-(P^. The integrals /„ (see 
appendix) contain two parts, the infinite T = contribution and the finite T > contri- 
bution. The infinite T = part requires a renormalisation procedure. In a strict sense, the 
starting Lagrangian (^ is not renormalizable because of the logarithmic interaction po- 
tential (|l]). But a possible regular isat ion of the infinite integrals would be the introduction 
of a cut-off parameter A. This regularisation is often used in effective non-renormalizable 
theories, where the cut-off determines the scale up to which the effective description is 
valid. The standard example is the Nambu and Jona-Lasinio (NJL) model. 

But such a kind of additional parameter violates the scaling properties of the effective 
theory, which we have chosen to simulate QCD, because of the appearance of a new scale. 
In order to keep the wanted scaling behaviour one has to give the cut-off parameter also a 
conformal weight of one unit by some mechanism, as for example done for the NJL model 
in Si- 

We simplify the procedure here and cancel all infinite T = contributions by hand. 
This means that we normalize the theory at tree level. We will determine the parameters 
C and Xp in (0) such that the EP at T = reproduces the bag constant B and the glueball 
mass mx'- 

B = V,^{X = 0) - V,^{X = Xo) , (33) 

ml = Ml{T = 0) = V::^{X = X,) . (34) 

As we will see later the determination of the EP at the origin A? = is not possible because 
of infinities in the effective coupling constant in the expansion of the logarithm. Thus a 
renormalisation, i.e. a determination of the parameters via ( p3| ) is not possible. 

By applying (^) and (0), we fix the depth and the curvature at the minimum of 
the EP at T = 0. We are dominantly interested in the thermal behaviour and therefore 
in the thermal excitation around this minimum. But the characteristic of this behaviour 
is already mainly fixed by these two parameters (depth and curvature). Therefore we 
expect no drastic change in the result of the thermal properties by changing the level of 
normalisation in the theory. 

For the value of the bag constant we choose as two limiting cases B^^^ = 140 MeV and 
^ 240 MeV which bound the common range of values appearing in the literature. 
We want to emphasize that the larger value is more realistic as newer results from hadron 
spectroscopy show ||2^, compared to the original work of the MIT-bag model where 



the lower B was used. Also, only the larger value is compatible with the gluon condensate 



at T = as extracted from QCD sum rules 



For the glueball mass we choose mx = 1-6 GeV, but we will also investigate other 
values. That value is motivated by lattice calculations |^, and also the experimental 
search favours candidates in the mass region of 1.5-1.8 GeV pB|. 
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For the parameters of the chiral sector (c, rriq) we choose standard values which repro- 
duce the physical pion mass and, for the case of SU(3)b, also the K and t] masses in a 
reasonable way. We set f^, = 93 MeV. 

We now apply to the partition functions (P7D-(P^ and neglect as discussed 

-"int 



the T = contributions in the integrals /„. Up to first order in the expansion of we 



are allowed to set 29 



X = X,i. (35) 
For the case of a pure dilaton field theory we get for the GEP 

v^^ix) = v(x) + If (fi2, T) + liv"(x) - ^]2)/o^(^]^ r) + ^v""(x)iUn^ rf , (36) 

and the variational equation for the determination of VL reads 

\{v"{x) - n') + \v""(x)iUn\ t) = o . (37) 

The result is plotted in Fig. 2 for various temperatures and for the two different bag 
constants. In that figure the singular region at A" = has been truncated. We see from 
eq. (0) that the four-point coupling of the .Y-field, which is proportional to V""{X) = 
24 C [in [X/X) + 11/6], develops a singularity at X = 0. The expansion of the tree level 
potential breaks down at the origin. According to the interpretation of the model, the 
phase X = would correspond to a vanishing gluon condensate, and because of eq. (^) 
also to a vanishing chiral condensate. This means that all of the non-perturbative physics 
has gone and we are in the pure perturbative region. An interpretation of this phase as 
the deconfinement phase as done in is therefore natural. But this would mean that 
the effective degrees of freedom of our model, mesons and glueballs, have been dissolved. 
Therefore a breakdown of the model at A* = is only to be expected. 

However, we can still use Fig. 2 to extract the temperature scale, where thermal ex- 
citations become important, even if we can't explicitly follow the phase transition to 
X = 0. This temperature scale is apparently determined by the chosen bag constant. For 
B^/^ = 140 MeV we begin to see a strong shift of the minimum at temperatures of about 
T = 250 MeV, and above T = 300 MeV the minimum at A" 7^ is lost, while for the large 
bag constant B^^'^ = 240 MeV the GEP doesn't show any visible shift up to temperatures 
of order T = 350 MeV, and only above T = 450 MeV the minimum ai X ^ disappears. 

We now proceed to the GEP of the full Lagrangian (^, where we choose a SU(A^) 
flavour symmetric mass matrix Aig. The result for the GEP is 

vS'^''\x) = v''(x) + iUn',T) + h{vy(x)-n']i^in\T)+ lv""(x)i^in\TY 



+ iN' - (M^iX), T) + ^ ^^^^^ - (^^ j U {M^{X), T) 

+ ^-^^^^^i^iM^A-^), T)I^{n\ T) , (3^ 
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with the condition for 



1 (N'^ 



^li{M^{X),T)mn\T) = 0. (39) 



The result is plotted for = 2 in Fig. 3. Qualitatively the GEP behaves similary to the 
pure dilaton theory, Fig. 2. The dominant contribution to the shape of the GEP (^) as 
a function of X comes from the one loop term IJ'[fl'^{X), T] of the glueball field. (This 
does not imply the same for the absolute values of the different terms in ( |38D , because the 
overall normalisation of the curve in Fig.'s 2, 3 was chosen by hand to facilate comparison 
at different temperatures.) This means that the glueball dynamics itself dominates the 
position of the minimum of the GEP. The Goldstone bosons and their coupling to the 
dilaton field play only a minor role. 

The leading singularity at the origin, however, is influenced strongly by the four-point 
coupling of the mesons and the coupling between glueballs and mesons which diverge as 
1/X as X 0. Therefore the two last terms in ( |38D dominate the singular behaviour 
and change the sign of the pole (now positive) relative to the pure glueball case. These 
explains the steep rise of the EP curves near the origin at high temperatures. 

As a final example we investigate the case of broken flavour SU(3)b. The result for the 
GEP is 



v: 



SU(3)f 



cS 



Vix) + jf (^]^ T) + l[ivy(x) - n']i^in\ t) + \v""{x)i^{n\ Tf 

M^{Ml(X), T) + AI^{MWX), T) + I^{Ml(X), T) 
1 /^o 



24^ V X 



<dmll^{Ml{X), Tf - {ml + Acm,)I^ {Ml{X), Tf 

- Qmll^{Ml(X), T)I^{M'^(X), T) 
+ lQmll^{Ml(X), T)I^{M'^(X), T) 



'2X Xn 



U {Sl\ T) SmiU (M^iX), T) + 4m^/o^ (M^(A'), T) 



+ m;/o^(M^(^),T) 



with 



+ 



{X) - 



+ -^v""(x)i^{n\ T) 



(40) 



2,X Xq 



?>mm (M^^(^), T) + Amy^ (M^(A'), T) + m'U {M'{X), T) = .(41) 



Here we only have a result for the large value of the bag constant, shown in Fig. 4. The 
reason is that for the potential V'^{X) (^ the negative curvature due to the g^X^-teim 
is so strong that by varying C and Xp the bag constant can never be reduced below the 
value S^/^ = 220 MeV. The general behaviour of the curves in Fig. 4 is similar to the 
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former cases. The onset of the thermal shift of the minimum occurs at somewhat lower 
temperatures; this can be seen better below, when we consider the gluon condensate. 



3. CONDENSATES 

3.1. The Gluon Condensate 

The minimum of the GEP determines the vacuum expectation value {X). This minimum 



where we use the approximation (O") ~ (O)", because the EP allows only the determina- 
tion of {X). The result is shown in Fig. 5 for different parameter sets. 

The crucial point is that the gluon condensate is stable up to temperatures of order 200 
MeV for both bag constants. We already mentioned that only the higher bag constant is 
realistic, where the decrease of the gluon condensate sets in even later. However, we must 
emphasize that our results should be taken only on a very qualitative level above T = 150- 
200 MeV. The reason is that our starting point, the non-linear sigma model, represents 
only the lowest order of chiral perturbation theory, and the neglecting of higher order 
gradient terms limits its quantitative applicability to temperatures below the pion mass. 
Thus the relevant change in the gluon condensate happens only at temperatures where 
the chiral sector of the model has already broken down. However, the critical temperature 
for the melting of the gluon condensate is dominated by the glueball dynamics, as can be 
seen from a comparison in Fig. 5 of the pure glueball scenario with the model including 
all eight pseudoscalar mesons. For this reason we believe that our statements about the 
behaviour of the gluon condensate at temperatures above the limit of validity of chiral 
perturbation theory are at least qualitatively correct, unless the chiral phase transition 
implies also deconfinement for the glueballs (for which our effective model does not provide 
any mechanism). 

We see in Fig. 5 that the value for Tc is dominated by the value of the bag constant. 
A larger bag constant results in a higher Tc. At the high-temperature end of the drop 
in the gluon condensate the curves seem to level off; thus the gluon condensate doesn't 
completely vanish at Tc, but approaches zero only in the limit T — > cxd. The surviving 
condensate just above the steep drop seems to be largest for the case of SU(3)b with a 
large strange quark mass. This ties in with the observation made in the introduction that 
the gluon condensate is strongly correlated to the scale anomaly. The anomaly, a quantum 
effect, is not expected to vanish at some fixed temperature, in contrast to spontaneously 
broken symmetries |T^, |3^ . What we apparently see here is the fading of the condensate 
relative to a new scale brought in by the temperature. 

To be more specific, let us for simplicity look at the pure glueball case. The scale of 
the anomaly is given by the gluon condensate at T = or by the bag constant, which are 
related in our pure glueball model by 




(42) 



B 



(43) 



4^ 2g 
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We now read off the critical temperatures for tlie two investigated bag constants, obtaining 
% = 290 MeV for = uq UeV and T, = 470 MeV for B^/^ = 240 MeV. We tlius have 



Tc ^ 2 X B'/' = V2 {(Mc-^C^nT^o ■ (44) 

The result = 0{{G^) 

T=o) is no surprise, since there is no other relevant scale. (The 
glueball mass is of minor influence, as will be shown below.) However, it confirms the 
interpretation, discussed above, that at Tc the temperature becomes the dominating scale, 
leading to a fading of the gluon condensate. But no complete vanishing is observed, 
meaning the anomaly stays on at high temperatures. Therefore the non-perturbative 
physics incorporated in the gluon condensate doesn't vanishes at but stays on to very 
high temperatures. A continuation of non-perturbative physics beyond the of the 
Wilson loop is also seen in other sectors, e.g. instanton effects ||31| or chromo-magnetic 
correlations p^ . 



If we look at the influence of the chiral dynamics on the behaviour of the gluon con- 
densate, we see only minor effects. The general behaviour is that as more mesons are 
incorporated as more is lowered, but the corresponding shift is very small. A larger 
effect is seen when we include flavour symmetry breaking in the SU(3) case by turning on 
the strange quark mass and bringing the kaons and rj to their physical mass. This effect 
is at first surprising, because naively the heavier mesons shouldn't be excited as easily as 
in the case of exact SU(3) symmetry with (nearly) massless quarks. The reason for the 
observed effect is the scaling behaviour of the meson masses (^I[) , as we will now discuss. 

We determine the gluon condensate via the minimum of the EP. We vary X and 
determine the energy density. Important are, however, not the absolute contributions, but 
the variation with X, because the global normalisation of the EP is arbitrary. All meson 



masses scale with \J X / Xq. Therefore in the thermal weights we encounter expressions like 

g-mV^/T ^ (45) 

which vary more strongly with X the larger the rest mass m is. Thus the location of the 
minimum is more sensitive to the heavier chiral mesons. 

This, of course, raises the question of the influence of all the neglected other heavy 
hadrons. They surely are important for temperatures of order 200 MeV and higher. 
Their influence would probably be a shift of to lower values, as indicated above. But 
the incorporation of these other states into the model is problematic because for these 
non-Goldstone particles the scaling properties with X are not known. This limits the 
quantitative power of predictions of our model in the high temperature region. 

We summarize by stating that the gluon condensate is very stable up to tempera- 
tures above 200 MeV. This is in agreement with other results for the temperature (in- 



) dependence of the gluon condensate within effective low-energy models |p!0| , |TT[. The only 
exception is the work by Bernard et al. 0, but this is probably due to the very low bag 
constant taken in their approach. 
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We now vary the glueball mass at T = in order to study the influence of this second 
empirical value on our results. The gluon condensate as a function of T is shown in Fig. 6 
for three values of mx- We see that a lower glueball mass drives Tc to lower values, 
since lighter glueballs are more easily excited at a given temperature, thus melting the 
condensate earlier. Quantitatively, however, the effect is small: a variation of the glueball 
mass by 1 GeV (from 1 GeV to 2 GeV) shifts Tc only by about 50 MeV. Thus the dominant 
scale for Tc is indeed the bag constant vis. the gluon condensate at T = 0. 

3.2. The Quark Condensate 

We finally investigate the influence of the coupling between dilaton and chiral fields on 
the chiral condensate. In QCD the chiral condensate is obtained by taking the derivative 
of the partition function with respect to the current quark masses. 



We replace the QCD generating functional by the effective partition function (P7|)-(PUD. 
We can use eq. (p6|) and rewrite the derivative with respect to into a derivative with 



respect to m^. This is equivalent to using the Gell-Mann-Oakes-Renner relation 

2/>2 = _(^^ + md){0\uu + dd\0) + 0{ml) , (47) 

which connects the QCD parameters to the empirical parameters /tt, and {qq) at 
zero temperature. If we apply djdmT^ to (p7D-(p9|), we need to account for the implicit 
dependence of on via eq. (^, which results in the equation 

{y^)\Xc{) = . (48) 

Using this condition we get 



Therefore we can neglect this dependence. Similarly the variational equation ( P^ causes 
a dependence of on m^r. It is of the same order as (|49|), and therefore we neglect it, too. 
We thus get for the chiral condensate 

(gg)^^(^)(r) ^ 

(gg)(r = 0) Uo, 



^-^(^^o(M.W,T) 



2N^f^ \Xci 
jN' - l)X,' 



iUm'A^ci), T)ll{^\ T) + -2 I^{n\Xc{), T) 
I^{Ml{Xc{), T)l\{Ml{Xc{), T) 

(50) 



2mf^ \Xci 
+ ^'^2j^;j^3^^° ^^(^^ T)/T,(M^^(^.0, T) 
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If we set (Xci/Xq) = 1 and cancel all terms containing the dilaton mass parameter Q, 
we recover the well known result of chiral perturbation theory |]3^. The corrections due 



to the dilaton field are of two kinds. First we have the universal scaling by {Xci/Xq)^ as 
already shown in eq.(^). Second also the {Tt{U + f/^)) factor receives corrections because 
of the coupling between the dilaton and the chiral field. First there is a scaling of the 
terms coming from pure meson loops because of the scaling of the meson propagators 
with negative powers of (Xci/Xq) (plQl). Second there are contributions from glueball loops. 



These last terms are small as has already been shown in ||T^ . 

The quantitativly largest modification comes from the universal scaling, which goes 
with the third power of X^. The result is plotted in Fig. 7 for a flavour SU(2), where 
we set Xci = X and use the temperature dependence of X and fl as extracted from 
the GEP. The solid line corresponds to the result of chiral perturbation theory ||3^. The 



modifications are really small. For the realistic larger bag constant there is nearly no visible 
shift, because [X/Xq) ^ 1 in the whole temperature region where the chiral condensate 
is non-zero. Only the smaller bag constant leads to a downward shift of order 20 MeV of 
the chiral T^. In general we conclude that the scaling properties of QCD have only a very 
weak influence on the chiral dynamics. 



4. CONCLUSIONS 

We have started from an effective Lagrangian which implements the chiral symmetry and 
scaling aspects of QCD. This allows the investigation of two interesting questions, the 
thermal behaviour of the gluon condensate and the influence of the scaling properties on 
the well investigated chiral restoration phenomenon. 

We find that the gluon condensate is very stable up to temperatures of 200 MeV, where 
the chiral sector of the theory reaches its limit of validity. As a result the chiral dynamics 
is hardly changed at all up to the chiral phase transition. Thermal variations of the gluon 
condensate and of the chiral condensate occur at two quite different temperature scales. 
A conclusion from our work is that the gluon condensate does not drive the disappearance 
of the chiral condensate at high temperature as suggested as one possible scenario in ^. 

In our model the melting of the gluon condensate is dictated by the glueball dynamics 
itself. The influence of the chiral mesons is minor, but other hadron states neglected in 
our approach may further change the temperature scale of the gluon condensate. This 
thermal scale is fixed by the value of the condensate at zero temperature or, equivalently, 
the bag constant. In a pure glueball theory we find Tj'^'^'* ^ y/2{GG)]^*Q. 

It was suggested in that the expectation value {X) of the dilaton field, i.e. the gluon 
condensate, could be used as an order parameter for gluon deconfinement. Our results 
contradict this interpretation. We see a relatively stable gluon condensate, well beyond the 
point of evaporation of the quark condensate which, according to |^ , should be interpreted 
as quark deconfinement. Since we see no reason why one kind of colored particle should 
be deconfined earlier than another one, we suggest that the color deconfinement phase 
transition occurs before the gluon condensate vanishes. This implies that nonperturbative 
phenomena persist well into the deconfined phase. 

This is also seen in lattice QCD results, where one is able to extract via some extra- 
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polation the finite temperature behaviour of the gluon condensate ||35|. There one sees at 
the critical temperature of the Wilson loop no evidence for a rapid decrease of the 
gluon condensate. The authors of ||35| suggest that at half of the zero temperature 
condensate survives. All of these analysis suggest that the gluon condensate is not suited 
as an order parameter for deconfinement. 

This conclusion is also supported by the fact that the gluon condensate is strongly con- 
nected to the QCD scale anomaly. Quantum-field theoretical anomalies are not expected 
to vanish at some finite temperature or density. Our results for the gluon condensate seem 
to support this expectation, since in our calculation the condensate, after a steep decrease 
at some T^, levels off at a finite value, suggesting complete disappearance only in the limit 
T oo. 
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Appendix 

Here we present a list of integrals used in the text. We use the notation of |jl9|. The zero 
temperature integrals are defined by 



2N 



2E{p){2Txf 

These and their finite temperature analogues satisfy the following recursion relation: 

In this paper we have used the following integrals: 

J_i(M^ T) = /Tr°(M2) + /Ti(M^ T) 



/Ti(M2, T) 



d^p 1 



(51) 

(52) 

(53) 
(54) 



(27r)3E2(p) {e(^m - 1) 



1 f3e^^(f^ 
+ 



E{p) e/5^(P^-l 



1 



47r2 



dx 



x'^ie^ — 1)^ 



[(1 + x)e" - 1] 



/o(M^ T) 



f3M 

I^=\M^) + /if (M^, T) 

d^p 1 
(27r)3 2E{p) 



(55) 

(56) 
(57) 
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° ^ ' ^ i (27r)3 e/3^(^^ - 1 27r2 i - 1 



I^{M^ = 0, T) 



j'2 

12 

/i(M2, T) = lJ=\Ke) + lJ{M^,T) 
^ - J (27r)3 2 

OO 



27r2 

I3M 



li (M^ = 0, T) = 



90 
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Figure Caption 



Figure 1.: Feynman diagrams for tlie calculation of (0) in first order in the expan- 
sion of exp{i / d'^xC\l^}. The solid line corresponds to the glueball field X, the long 
dashed line to and the short dashed line to 9^$". 

Figure 2.: The Gaussian effective potential for the case of only dilaton fields for 
various temperatures T. The zero temperature glueball mass is m;^ = 1.6 GeV and 
B^/^ = 140 MeV in a) and B^^^ = 240 MeV in b), respectively. The singular region at 
X = has been truncated. 



Figure 3.: The Gaussian effective potential for the case of flavour SU(2) in the meson 
sector for various temperatures T. The zero temperature glueball mass is mx = 1-6 
GeV and B^/'^ = 140 MeV in a) and S^/^ = 240 MeV in b), respectively. The singular 
region at A" = has been truncated. 

Figure 4.: The Gaussian effective potential for the case of broken flavour SU(3) 
{ttIu = fUd < rris) in the meson sector, which represents the physical meson masses, 
for various temperatures T. The zero temperature glueball mass is mx = 1.6 GeV 
and B^/^ = 240 MeV. The singular region at A" = has been truncated. 



Figure 5.: The gluon condensate as a function of temperature, normalized to its zero 
temperature value. The solid line corresponds to the scenario of glueball fields only, 
the long dashed line to flavour SU(2), the long-short dashed line to flavour SU(3), 
and the short dashed line to broken flavour SU(3) (see Fig. 4). We show results for 
two bag constants and mx = 1.6 GeV. 

Figure 6.: The gluon condensate as a function of temperature, normalized to its 
zero temperature value. For the meson sector we use an SU(2) symmetric spectrum, 
and we set B^^'^ = 140 MeV. The solid line corresponds to mx = 2.0 GeV, the short 
dashed line to nix = 1.6 GeV and the long dashed line to m^f = 1-0 GeV. 



Figure 7.: The chiral condensate as a function of temperature, normalized to its 
zero temperature value. For the meson sector we use an SU(2) symmetric spectrum. 
The solid line corresponds to the result of chiral perturbation theory |Q . The dashed 
lines corresponds to the corrections due to the dilaton field. The short dashed line is 
for iS^/^ = 240 MeV and the long dashed line for B^/^ = 140 MeV. 
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